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Abstract: We investigate to what extent the recently measured value for a non- vanishing 
^ H direct CP asymmetry in D° — > K + K~ and D° — > 7r + 7r~ decays can be accommodated 

in the Standard Model (SM) or extensions with a constrained flavour sector, for instance 
from a sequential 4 th generation of quarks (4G). From the comparison with D° -> K~ 
branching ratios, we establish large U-spin symmetry (d 4-¥ s) breaking effects with large 
strong phases between different interfering amplitudes. On the basis of conservative esti- 
mates on amplitude ratios — which are supported by an analysis of the breaking ofacffu 
symmetry in non-leptonic B° decays — we find that, in the SM, direct CP asymmetries in 
the 7r + 7r _ or K + K~ modes (or in their difference) of the order of several per mille are still 
plausible. Due to the constraints on the new CP phases in the 4G model, only moderate 
effects compared to the SM estimates are possible. We suggest CP studies at LHCb as 
well as at (Super)B-factories of several distinctive modes, such as D + — > K~W°7r + , (fyn + and 
7r + , <f)ir + (K + ) etc., which should shed more light on the short- and long-distance 
issues underlying CP violation in non-leptonic D-meson decays. 
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1 Motivation 

CP-violation studies in the charm system are important for a variety of reasons. First, 
in the Standard Model (SM), only small asymmetries are expected, see e.g. [1-4] and 
references therein. On the other hand, in many models of new physics (NP), the top 
quark is sensitive to non-standard effects, and there are flavour- changing interactions which 
couple the charm and the top quark. Flavour models based on generic ideas of warped 
extra dimensions represent a very interesting example of this type; especially so as large 
0(1) CP-odd phases may accompany the new interactions [5]. This is also the case in many 
models with an extended Higgs sector [6]. Even in a simple extension of the SM with a 
fourth generation of fermions, significant differences in CP asymmetries in the charm sector 
may arise as compared to the SM [7-9]. Furthermore, the charm-quark mass (~1.2 GeV) 
is not that heavy; consequently large (CP-even) rescattering phases can be present in D- 
meson decays which, together with CP-odd phases from the SM and/or NP sector, can lead 
to sizeable direct CP asymmetries. CP studies via D°-D° mixing are also highly motivated 
as the D° is a unique bound state of charge 2/3 quarks and is thus sensitive to NP affecting 
the top quark. 

In this article, we will be mainly concerned with CP violation in non-leptonic D° 
decays. The time-integrated CP asymmetry Acp(f) for a given final state / is defined as 



Acp(/) 



F(D° -> /) - T(D -> /) 

r(D0_»/) + r(5°-»/) 



(i.i) 
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which (to first approximation) may be decomposed as (see e.g. [10]) 



A CP (f) = A d £ + A^ + A™. (1.2) 

Here, A^,(f) is the direct CP asymmetry in the decay D° — > /, A^p is the CP asymmetry 
from D°-D° mixing, and A™p stems from the interference of mixing and decay. Recent 
results from the LHCb experiment [11] on CP asymmetries in D° decays, 

A-Aqp = A$p(K + K-) - ^(vr+vr-) = - (0.82 ± 0.21 ± 0.11) % , (1.3) 

indicate a 3.5a deviation from 0, with a large amount of experimental systematics cancelling 
in the considered difference of decay modes. We remind the reader that, in the SU(3)p 
limit for light quarks [10], the direct CP asymmetries in the individual channels are equal 
in magnitude but opposite in sign, Aq P (K + K~) ~ — Aqp(it + it~). Assuming that the 
effects of Aqp drop out in the difference, and that the effect of A^p is small, this points 
to a relatively large value for |.Aqp| for both decay modes. Not surprisingly, the LHCb 
result has renewed the interest in CP violation in charm physics, and a number of articles, 
addressing the interpretation of the data within the SM or NP, have appeared since (see 
e.g. [12-23]). At about the same time, the CDF collaboration has reported their result 
from an independent measurement [24] , which also provides numbers for the individual CP 
asymmetries in the (K + K~) and (ir + ir~) channel, 

A C p(K + K-) = (-0.24 ± 0.22 ± 0.09)% , 
A C p(vr+7r") = (+0.22 ± 0.24 ± 0.11)% , (1.4) 

which is consistent with the above LHCb result. Both measurements dominate the world 
average provided by the HFAG collaboration [25] , 

AAqp = (-0.645 ± 0.180)% . (1.5) 

Although a direct CP asymmetry of the order of a few per-mille seems to be a small 
number, the SM expectations for Ai^p in D decays are usually even smaller, because 
the only source of observable CP violation comes from the interference with a sub-leading 
amplitude which involves a relative suppression factor involving 4 powers of the Cabibbo 
angle. A precise quantitative estimate, however, is difficult because standard approximation 
methods known from i?-meson or kaon physics do not seem to work sufficiently well in the 
charm sector, since the charm-quark mass is neither much smaller nor much larger than 
the typical hadronic scales. Without strong theoretical prejudices about non-perturbative 
hadronic dynamics, we are therefore constrained to semi-quantitative analyses. The related 
questions that we are going to explore in this paper are: 

• What can be learned from U-spin relations between down- and strange-quarks, and 
how large is the effect of U-spin breaking? 

• Do we have phenomenological evidence for large strong phases in D° — > P + P~ decays, 
and how does this compare with the situation in non-leptonic _B°-meson decays? 
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• How large do sub-leading amplitudes in the SM have to be in order to explain the 
observed AAg$? 

• To what extent can NP models with constrained flavour sector, like for instance a 
model with 4 fermion generations, enhance the theoretical predictions for A^Qp as 
compared to the SM? 

• Can we identify other D-decay modes which should exhibit similar patterns of CP 
violation? 

2 [/-Spin in Non-Leptonic D° — >■ P + P Decays 

Assuming SU(3) flavour symmetry for the strong interactions of light quarks, the con- 
tributions to the amplitudes for non-leptonic D-meson decays can be related. For our 
purposes, it will be sufficient to focus on a sub-group, U-spin symmetry, which acts on 
U-spin doublets (d, s) and which have been frequently used to analyse weak non-leptonic 
meson decays, see for instance [26-29]. 

Starting point is the weak effective Hamiltonian for c —> uqq' transitions (q, q' = u, d, s) 
which can be decomposed into a U = and U = 1 part, as follows 

H e g(c — > usd) 

H eS (c -> uqq) 

H c s(c —t- uds) 

Here Hu = \ only contains the current-current operators from tree-level TU-boson exchange, 
while Hjj = q also receives contributions from strong and electro-weak penguin operators. 
The important property to notice is that the U-spin singlet or triplet terms are multiplied by 
a single combination of Cabibbo-Kobayashi-Maskawa (CKM) elements. In the exact U-spin 
limit, the amplitudes for the related Z)°-meson decays into two light charged pseudoscalar 
mesons 1 are thus described by only two independent complex amplitudes, which we denote 
by A U=Q and B u=1 , 

2 V^Vud Bu=i , 

(Ad + A s ) Au=o + (Ad - \ s ) Bu = \ , 
(Ad + A s ) Au=o — (Ad — As) Bu=i , 
2V: d V us B u=1 , (2.2) 

where we have defined Ad = V* d V uc i and A s = V* S V US . From the Wolfenstein expansion of 
the CKM elements in powers of A = sin#c ~ 0(0.2), we infer that the decay D° — > K~tt + 

*We are restricting ourselves to charged mesons (with U = 1/2) here, because in this way we only 
generate U — and U = 1 final states which require a minimal set of independent amplitudes, once we 
allow for arbitrary U-spin breaking, see below. 



fv; d v ud -v*v us 



V V2 

(ou^Ct +1) 



i4=r 0) + (KdKd + v*v us ) h u=0 , 



(2.1) 



U-spin limit: A[D° -)• K~ir + ] = 

A[D° TT+TT-] = 

A[D° -> K + K-} = 
AID -»• K + tt-] = 
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BR[Z?° - 


-> K 


-7T+] 


BR[L>° 


->■ 7T" 


NT] 


BR[D° - 




K~] 


BR[L>° - 


» K" 


I-7T-] 


BR[D° - 


» K- 


-7T+] 



is Cabibbo allowed (CA), the decays Z)° — > 7r + 7r~, K + K~ are singly-Cabibbo-suppressed 
(SCS) with (Ad — A s ) ~ O(A), while the decay D° — > K + ir~ is double-Cabibbo-suppressed 
(DCS). Furthermore, the contribution of the U = amplitude to the decays D° — > 7r + 7r~ 
and D° — > K + K~ is suppressed by another 4 powers of A, since A^ + A s ~ 0(X 5 ), and 
therefore the decay rates should be equal, while the CP asymmetries from the interference 
of Ajj = q and Bjj = i should be tiny, as already mentioned above. However, the measured 
branching ratios (BRs) for the CA, SCS, DCS modes (experimental numbers are taken 
from [25] or [49], see also references to the original experiments therein) 

(3.949 ± 0.023 ± 0.040 ± 0.025)% , 
(0.1425 ± 0.0019 ± 0.0018 ± 0.0014)% , 
(0.3941 ± 0.0038 ± 0.0050 ± 0.0024)% , 

(0.331 ±0.008)%, (2.3) 

do not follow the pattern expected from U-spin symmetry. Actually, correcting for phase- 
space effects and CKM factors, the following observables 

obS3S B?pr^7F^ A ±0.03, (2.6) 

deviate significantly from unity. In principle, this can be due to: U-spin violation origi- 
nating from m s ^ rrid in the SM, and/or a drastic enhancement of the hadronic matrix 
elements of penguin operators in Ajj = q either from long-distance dynamics in the SM or 
short-distance effects from sizeable NP. In the latter case, however, it can easily be seen 
that the enhancement of order A -4 ~ O(600) necessary to reproduce the pattern of BRs, 
at the same time clashes with the (still small) values for the direct CP asymmetries, which 
in such a case would be enhanced by the same factor, unless the strong phases of Ajj=q 
and -B[/=i were fine-tuned. Actually, we will find below that the simultaneous fit to obsi 
and obs2 clearly fixes the strong-phase differences of the relevant interfering amplitudes in 
a scenario with broken U-spin to be large. 

In the following, we are thus allowing for generic U-spin violation in the above decay 
amplitudes. In order to explain the measured BRs, we would have to accept U-spin violating 
effects as large as 50-60% on the amplitude level. In the naive factorization approach, such 
factors can largely be explained by the difference in decay constants (fx/ fn) and hadronic 
form factors {F{D — > K)/F(D — > tt)) each of which amounts to O(20%) corrections to 
the U-spin limit (recent lattice results 2 on these quantities can, for instance, be found in 
[30, 31]). As already mentioned, the observed BRs also require a significant strong-phase 



2 We thank Ruth Van de Water for discussions pertaining to these. 
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difference which, together with the remaining amount of U-spin violation in the magnitude 
of the decay amplitudes, points towards essential non-factorizable long-distance effects in 
the hadronic matrix elements. Again, this comes as no surprise, as similar conclusions 
have been drawn from the analysis of non-leptonic 5-meson decays (see e.g. [32-42] for an 
incomplete list of references), although the importance of non-factorizable contributions 
is appreciably less dramatic in .B-meson decays due to the fact that the expansion 
there is more efficient than a l/m c expansion in D-meson decays. 

2.1 Including U-Spin Breaking 

A non-zero strange-quark mass in the strong-interaction Hamiltonian gives rise to a new 
U-spin triplet operator that enters the hadronic interactions. In particular, this can turn 
U = operators in the weak effective Hamiltonian into U = 1 amplitudes (denoted as 
ABfj =1 ) in the hadronic matrix elements and vice versa (Hu=\ contributes to AAu=q). 
The decomposition of the relevant decay amplitudes including (first-order) U-spin breaking 
then can be written as (see also [18, 23]) 



A[D° -> K-ir + ] = 2 V* s V ud (B U=1 - AB' U=1 ) = 2 V c * s V ud B u=l 



1- tie* * 



A [ D o ^ Tr+^-j = (Ad + Xs ) (A u=0 + AB U=1 ) + (Ad - X s ) (B u=1 + AA 



u=o) 



(Ad + A s ) re^ + ne** 1 +(Ad-A a ) l + r e 



= B V -_ 

A[D° -> K+K~] = (Ad + A.) (A u=0 - AB U=1 ) - (X d - X s ) (B u=1 - AA 



J <t>0 



Bu, 



U=0) 



(Xd + X s ) (re^-ne^ ) - (X d - X s ) l-r e^° 



X[D° -> K+7T-] = 2 v: d V us (B u= i + AB' U=1 ) = 2 v: d V us B u=l [l + r[ e^] . (2.7) 

In the second equations of this parametrization, we have introduced various amplitude 
ratios and strong phases with respect to the originally leading By = \ amplitude. 

It is important to note that in the SM — or for that matter in many NP models - 
direct CP asymmetries in D° — > K^tt^ cannot arise due to the absence of penguin contri- 
butions in these decay modes. However, D and D° both can decay to the same final state 
K^tt^ which leads to interference between D°-D° mixing and the decay processes. Such 
a contribution though is expected to be extremely small (see below). On the other hand, 
if there are charged scalars with non-trivial flavour couplings and CP-odd phases (like, for 
instance in general models with extended Higgs sector), then direct CP asymmetries may 
arise in D° — > K^ir^, too. The largish BR of about 3.8% for D° — > K~n + thus provides a 
valuable opportunity to search for a non-standard CP phase. LHCb and Super-B factories 
producing about 10 9 charm mesons should have a 5a reach of searching for CP asymme- 
tries at the level of about 0.2% in this mode (similarly, for D° — > K + ir~ with a BR of 
about 1.5 • 10 -4 one can search for 0(1%) CP asymmetries). Note also that this simple 
decay mode is of crucial importance in extracting the CKM angle 7 from the B —> DK 
decays in the ADS [43] analysis, where traditionally one assumes no CP asymmetry in the 
subsequent D-decays. A confirmation of this assumption by direct experimental searches 
is therefore very valuable in any case. 
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Figure 1. Fit result for the amplitude parameters ro and cos fa (upper row) and r[ and cosfa 1 
(lower row), determining the amount of U-spin breaking in D° — > P + P~ BRs. The generic points 
shown in light blue are consistent with the experimental constraints at the 2cr-level and obey x 2 < 6. 
The black points denote a subset of points where the strong phase differences between Ajj=o and 
AA[/ = oi a s well as between Bjj = i and ABjj = i are assumed to be equal within a few percent, 
(4> - fa) = {0, tt} and fa = {0, w}. 



2.2 Phenomenological Constraints on Amplitude Ratios 

In the following, we perform a fit to the 8 amplitude parameters (r, ro, r\, r[ and <p, cfio, fa, fa) 
to the 3 experimentally measured ratios of BRs in (2.4-2.6) and the difference of CP asym- 
metries in (1.5). We also take into account the fitted value [25] of the strong phase difference 
between the D° — > K^tt^ decays, Ai^ = 22.4°^ 1 ^'q . To this end, we generated random 
points which (i) lie within the 2a ranges for each experimental observable, and (ii) yield a 
total x 2 -value of less than 6. We have allowed arbitrary strong phases while the magnitudes 
of the various amplitude ratios are considered in a conservative range \rx\ < 8. The result 
is illustrated in figs. 1,2. 

From fig. 1, we observe that the parameters describing the U-spin breaking effects from 
the amplitudes A^4o and AB[ are rather tightly constrained, with 

r ~0.52, cos0 o ^-O.64, ri~0.19, cos<£i>0.18. 

The individual values for the magnitudes (r, r±) and phases (fa fa) related to the Cabibbo- 
suppressed operators are not fixed from the fit. A useful parameter to study the effect on 
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Figure 2. Fit result for difference of CP asymmetries AA^p (top; horizontal lines indicate the 
la and 2a experimental constraints), and the ratio of the sum and difference of CP asymmetries 
(center) in D a — > K + K~ and D° — > tt + -k~ decays as a function of f = \A 2 /2 + rf/2. The 
lower plot shows the correlation between the two asymmetries acp3 = A^p(D° — » K + K~) and 
acp2 = Aqp(D° —> 7r + 7r _ ), (the dashed line indicates the naive U-spin limit). Plot conventions as 
in fig. 1. 

the direct CP asymmetries is given by the average 

f = ^2/2 + rf/2. (2.8) 

This is illustrated in fig. 2, where we have plotted the difference of CP asymmetries AAq^, 
as well as their sum divided by the difference as a function of f. The following interesting 
observations can be made: 
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• The effect of U-spin breaking from r[ (entering the BRs for D° — > K t^) is relatively 
small as compared to the effect from ro (entering the BRs for D — > K + K~ and 
D° — > 7r + 7r _ ). Again, this can be partly understood in the framework of naive 
factorization, where in the former case the U-spin breaking in the D —> K(ir) form 
factors tends to compensate the effect in the decay constants fxi^), whereas in the 
latter case, the two effects tend to add up. 

• There is a clear correlation between AA^'p and the minimal value for the effective 
amplitude parameter f. To reproduce the experimental 2<r-range for AA^p, one 
has to require f > 2. It remains to be seen whether such large values are consistent 
with our theoretical expectations, given our restricted knowledge on non-perturbative 
strong dynamics in time-like hadronic processes. 

• The U-spin relation, T,A^ = A^{D° -> K+K~) + A^(D° vr+vr-) = 0, re- 
ceives corrections of the same order as AA^p itself, or even larger (although the 
case Sj4^ p = is not excluded). Consequently, the correlation between Aflp(D — > 
K + K~) and Aq^,{D q — > tt + tt~) can be quite different compared to the naive U-spin 
limit. 

• The qualitative results for the BRs and CP asymmetries as a function of the amplitude 
ratios is well represented by a subset of parameters (indicated by the black points 
in figs. 1,2) where the strong phase differences between the amplitudes Au = q and 
AAjj = o (as well as between Bjj = i and ABjj = i) are set approximately to zero or tt. 

2.3 Simplified Analysis of U-Spin Breaking 

As we have seen above, the strong phase differences between individual U = 1 (or U = 0) 
amplitudes does not play an essential role for the analysis of direct CP asymmetries in D° 
decays (they simply reflect the redundancy in the effectively 3-dimensional space of the 
constrained parameter space in the fit). To simplify the further analysis, we will therefore 
set 

0i = (O,tt), <p= (fo^o + ir). (2.9) 

For the description of the BRs, we can also safely neglect the terms proportional to 
(Xd + A s ) ~ C(A 5 ). With this approximation, the parameters ro and 4>q can be extracted 
from the observed ratios of BRs, obsi_3. The constraint from obsi alone translates into 
the inequalities, 

2 2 

0.28 ~ 1 - -== < r < 1 + . ~ 3.51 , 

yobsi + 1 yobsi — 1 

obsi — 1 , 

cos (f> < — ^ -0.53 . (2.10 

obsi + 1 

Notice that these bounds already constrain the possible values for obs2 and obs3. Including 
the average of obs2 and obs3, the central values of ro and 4>o take values 

r ~0.54, cos 0o ~ -0.69, (2.11) 
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which are consistent with the previous fit within the uncertainties. 3 

Concerning the direct CP asymmetries, the constraints from obsi_3 alone cannot 
restrict the Cabibbo-suppressed amplitudes Ajj=o or AB[/=l which only provide a tiny 
(< 1%) correction to the BRs. Within our approximations, and assuming that all ampli- 
tude ratios are of 0(1), we obtain 

A^{K+K-) - A^(n+n-) * • £ + ^ ( ~ r ^ + ^ ^ ~ ^ , 

CPV 7 CPV 7 |A d -A s | 2 l-2r2 C os2</» + r-4 



together with 

A^[D° -»• ET+iT-] + -> vr+vr-] ^ obsi - 1 



(2.12) 

~-53%. (2.13) 



-> K+K-] - A^[D° -> vr+vr-] " obsi + 1 

The latter equation shows again that — as soon as U-spin violating effects are included 
- the direct CP asymmetries in D° — > it + tt and D° — > K + K~ are no longer equal 
in magnitude. The particular value in (2.13) reproduces the central value of the trend 
apparent in fig. 2 (right); allowing for arbitrary strong phases, the effect can become even 
larger. The former eq. (2.12) states that the measured difference of direct CP asymmetries 
(considering the leading term in the Wolfenstein expansion) is a product of 3 terms, 

-a5 + a:" 



/weak = 4 Im 



2A 2 A 4 r ? «0.11% 



.Ad — A s 

1 + rg 
1 — 2vq cos 2(j) Q + Tq 

/strong = —r sin (j) + r ri sin (j) cos <j>i . ( 2 -14) 

In order to achieve the central value of the experimentally observed number, /strong should 
be about —5.3, which requires accordingly large values for r and r%. More moderate values 
for r, r\ are possible if we consider the la range, only (as we have seen before). As the 
direct CP asymmetries scale linearly with both the amplitude ratios r and n, there is 
no preference for whether the original amplitude Ajj = q or the effect of U-spin violation 
through ABu=i should be dominating the required value of /strong, although the effect of 
n is somewhat diminished by the factor ro < 1. 

To illustrate the result, we consider the average of amplitude parameters f together 
with 

r sin0 r 

In fig. 3, we plot the result for /strong together with the experimentally preferred range for 
two values f = {2, 5} as a function of t. As before, we can see that with r = 5 (dashed line) 
and for negative values of t the result falls comfortably into the experimental la range, 
whereas for r = 2 (solid line), the result is marginally consistent with the measurement at 
the 2a level. 



The result of this procedure is similar to a recent analysis of 5(7(3) flavour-symmetry breaking in non- 
leptonic D decays considered in [18], where the authors neglected what would correspond to our parameter 
r'i which has been motivated by dropping higher SU (3) representations in the effective Hamiltonian. 
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Figure 3. The function /strong, determining the size of AA*p in the SM, as a function of the 
amplitude ratio t for two different values of the averaged amplitude ratios, f = 2 (solid lines) and 
f — 5 (dashed lines), see eq. (2.15). The values necessary to recover the experimental la (2a) range 
in eq. (1.5) are shown in different shades of grey. 



2.4 A Note on Mixing-induced CP Asymmetries 

As described in eq. (1.2), apart from the direct CP asymmetry, the time- integrated CP 
asymmetry Acp(f) in the decay of D° to a given final state / also contains a part which is 
due to the interference of decay with mixing, the so-called mixing-induced CP asymmetry 
which is defined as 



(l + lfl/l 2 ) 2 !^ 



Here Rf = denotes the ratio of the amplitudes for a given decay and its CP conjugate, 



Rm = p gives the ratio of the parameters p, q in D°-D mixing, x = denotes the ratio 

of the mass difference and decay width in the D°-D° system, and the phase (p denotes the 
argument of the function Xj = 2 The present experimental bound on the D° — D° 
mixing parameters, i.e. the value of x and the phase and magnitude of R m are given by 
[25], 

x = (0.98t°J4) % , arg (?) = (-8.5^)° , \R m \ = 0.87^;^ . (2.17) 

For decay modes like D° — > K + K~, 7r + 7r~, which — in the SM — do not have 
interfering amplitudes with different weak and strong phases of comparable size, one has 
\Rf \ ~ 1 and avg(Rf) ~ 0. Hence, the phase <f> will be purely induced by D°-D° mixing. In 
such cases ^4^p can be treated as a universal quantity which does not depend on the final 
state mesons. In particular, its effect will drop out in the difference of CP asymmetries, 
AAqp, apart from a term correcting for the finite experimental cut on the proper decay 
time, relative to the life-time associated to the individual decay modes (see e.g. [44]). We 
have checked that the presence of U-spin violating effects in A 1 q^(K + K~) and Aqp(tt + tt~) 
have a negligible numerical effect of the order 10~ 5 on AAcp. 
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2.5 AAqp in the Presence of a Fourth Generation 

In the 4G extension of the SM, the presence of the heavy &'-quark gives an extra contribution 
to the penguin operators, and it also modifies the CKM elements, such that + A s + Ab = 
— Ay / 0. Both effects can be accounted for by including an additional U = operator 
proportional Ay in the effective Hamiltonian. The amplitudes in eq. (2.7) can then be 
generalized to the 4G case as follows, 

A [ D ° ^ = Bu=i (r e i4> + n e i<Pi) + ( Ad _ Xs ) (l + r e** 

A[D° -> K+K-] = B u=1 \(-A b ) (re** - r, <""') - (A,, - A.J ( 1 - /•„ r 

(2-18) 



-Ay (r 4 e^ 4 -rie^* 



(') (') 

where r\ and cj)\ parametrize the contribution of the new U = operator and its U = 1 
counterpart from U-spin violation in the corresponding hadronic matrix elements. As 
before, we assume that the corresponding strong phases are related to the other U = and 
U = 1 amplitudes. Then, the additional contribution to Ayl^p can again be factorized, 

AAqp ~ /a£7 • { /weak • /strong + /weak " /strong } 

= /At/ " /weak " /strong > (2-19) 



with 



fwelk = 4Im 



Ad — A. s 



2 sin 6>i4 sin 6*24 sin(#i4 — ^24) 



sm»i2 



/strong = ~ r 4 sin 4 + r r 4 sin O cos 4 . (2.20) 

Here, we have used the PDG-type parametrization for the 4G CKM matrix as in [45], and 
neglected higher-order terms in the Wolfenstein expansion. Notice that eq. (2.13) is still 
valid if our assumptions on the strong phases hold. 

From the analysis of various flavour observables in the kaon and i?-meson sector (see 
e.g. [45-48]), we know that the product of mixing angles 614 and 624 can be as large as 0(A 4 ), 
up-to C(A 3 ), and therefore, in principle, can lead to an enhancement of the CP asymmetries 
by up to one order of magnitude as compared to the SM (where sin #13 sin 623 ~ 0(A 5 )). 
However, it has also been found that for such large values of 4G mixing angles, the new CP 
phases <5j 4 are rather fine-tuned to values satisfying 614 ~ 824 [8, 45]. This implies that the 
4G CKM elements cannot lead to a substantial parametric enhancement (i.e. a suppression 
factor with less powers of the Wolfenstein parameter A) of CP asymmetries in D decays 
(contrary to a recent claim in [17], where — to our understanding — the phenomenological 
constraints on the 4G CP phases have not been taken into account properly). Still, the 
presence of the additional penguin operators can lead to a numerical enhancement which, 
however, is hard to quantify reliably without more detailed knowledge on the hadronic 
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Figure 4. The function /s t ^ ong , determining the size of AA^p in the 4G extension of the SM, 
as a function of the amplitude ratio £4 for two different values of the averaged amplitude ratios, 
f 4 = 1 (solid lines) and f 4 = 3 (dashed lines). The associated weak factor is set equal to the SM 
one, / w J? ak = /weak, and the SM amplitude ratios are set to r = n = 2. The left (right) plot 
refers to a choice of strong phases where the two SM contributions r and r\ combine constructively 
(destructively) to a value of —2.2 (—0.67). The values necessary to recover the experimental lcr 
and 2cr range in eq. (1.5) are shown in different shades of grey. 



amplitude parameters. To illustrate the result, we define again 



s , .j±W s r^in^cos^ 

V 2 r4 sin 04 

and show in fig. 4 the result for /s t f ong as a function of £4 for f± = {1,3}, while we fix 
f4G = /weak and r = ri = 2, for simplicity, and compare two different choices for the 
associated phases corresponding to constructive/destructive effects in the SM contribution 
(assuming an overall negative sign). One sees that A-A^p can receive some deviations from 
the SM with 3 generations, but without a quantitative theoretical estimate of the strong 
amplitudes, nothing more can be said. 

2.6 Generic NP Explanation of AA$ p 

The above results can easily be generalized to generic NP models, and the following state- 
ments can be made: 

• Both, additional NP operators in the weak Hamiltonian with U = or U = 1, can 
contribute to the direct CP asymmetries. 

• The new flavour coefficients can potentially enhance the corresponding factor /^ e ^ k as 
compared to the SM. As a consequence, the observed value for AAq p can be obtained 
with somewhat smaller values for the various hadronic amplitudes. 

• Even in the presence of NP, in order to simultaneously explain the ratios of D° decay 
BRs, while keeping the direct CP asymmetries sufficiently small, we have to require 
large U-spin violating effects. Therefore, in general, the direct CP asymmetries in 
7r + 7r _ and K + K~ will not be equal in magnitude. 

For more detailed NP analyses, see e.g. [21, 22]. 



- 12 - 



Mode 


BR 


Acp in % 


5(7 Reach 


D+ -> K S 7T + 


1.47 x 10" 2 


-0.52 ±0.14 [25] 


1 x 10 -3 


D s -»■ r/7T + 


3.94 x 1(T 2 


-6.1 ±3.0 [49] 
-5.5 ± 3.7 ± 1.2 [25] 


0.7 x 10" 3 




1.21 x 10~ 3 


6.6 ± 3.3 [49] 
6.53 ±2.46 [25] 


4 x 10~ 3 



Table 1. BRs and CP asymmetries in different charged D + and D s decay modes. The quoted 
(naive) 5<r reach for the sensitivity on Acp refers to 10 9 produced D + ,D S mesons at LHCb or 
future Super-B factories. 

2.7 Possibilities for Direct CP Searches in Charged Modes 

Given that the current focus is on direct CP violation in neutral D° decay modes, we 
note in passing that decays of charged meson, D + and D s , also offer many interesting 
and experimentally distinctive channels for such CP studies, with the advantage that D°- 

mixing is not involved. In table 1, we list a few 2-body decay modes of D + and D s 
mesons where a non-vanishing CP asymmetry has been reported with ^ 2a significance. 
With the foreseen statistics at LHCb and future Super-B factories, we find that interesting 
sensitivity to CP-violating effects in the SM or beyond, to the level of a few per mille, can 
be anticipated for these examples, too. 

The current measurement of the CP asymmetry in the CA decay D + — > Kstt + 
within the present experimental uncertainties — is consistent with the SM expectation from 
mixing (which is of order e# ~ A 3 ), while the direct CP asymmetry is suppressed 
by a factor of A 6 and therefore negligible [50, 51]. NP contributions to the direct CP 
asymmetries of the order of several per mille could thus alter the SM prediction significantly 
which could be tested in the future. In this context, it would also be useful to tag on the 
flavour eigenstate in D + — > K°ir + (or also D + — > K*°ir + ). In a similar fashion, in the 
SCS decay mode D s — > Kg-K + , a contribution from K°-K° mixing is unavoidable, but in 
this case also penguin effects can contribute to Aq P at the per-mille level, as shown in the 
preceding discussion. The current experimental central value is rather high, but only with 
a small significance of about 2a. Here too, the decay into the flavour state, D s — > K^°tt + 
deserves experimental attention. Finally, the CA decay D s — > r]'ir + is a pure tree with 
vanishing CP asymmetry in the SM. Evidently, if the current 2a hint for the latter decay 
mode gets verified, close to its current central value, that could be an important sign of 
NP. 

There are many other modes of D and D s which could be suitable for direct CP 
studies. Some examples of SCS modes involving penguins are: D + — > K + K*° , K* + K°; 
D + -> (/>7T + ,p 7r + ,7r + 7r°(7/); D s ->■ K + 4>(r]'), K°(K*°)it + and many more. There are also 
further CA modes (with no penguin contributions) such as D + — > A >0 (A > *°)7r + , D s — > 
(j)ir + (K + ) etc. As mentioned previously, the SCS modes with penguins and with typical 
BRs of a few per mille should allow definitive CP searches of O(0.5%), and for CA modes 
w/o penguins we could get to about 0.2%. 
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3 Insights from Non-Leptonic £>-Meson Decays 



In this section, we are looking for a correspondence between the considered D° — > P + P~ 
modes and suitable non-leptonic i3-meson decays. Obviously, exchanging the roles of up- 
and down-type quarks would lead us to consider a symmetry between u and c, in the 
following referred to as IV-spin. Of course, we expect such a symmetry to be badly broken, 
and rather than use the symmetry relations themselves — for the sake of this paper — we 
would like to study the amplitude ratios in a parametrization based on TV-spin (including 
its violation to first order), analogous to eq. (2.7). In particular, we would like to quantify 
the phenomenologically allowed values for the various amplitude ratios in such a case, which 
might give us some guide-line about what to expect for the corresponding D° — > P + P~ 
case. 

With this in mind, we can repeat our previous analysis for the decays 



D--K+, D + K~,K~ir + , D~D"< 



in a straight-forward manner. 4 The relevant effective Hamiltonian, describing the b 
transitions, can be decomposed as 



H eS (b 
H cS (b 
H cS (b 



scu) 
sqq) 
sue) 



(v cb v: s ) H ^r 1] 



v ub v:. 



VcbV* 



V2 



H 



(W 3 =o) 
w=i 



+ (VubVus + VcbVcs) H W -_ 



> sqq 

(3.1) 
(3.2) 
(3.3) 



Notice that, in this case, the CKM elements scale as 



A c = V cb v: s ~ AX Z , 

K = v ub v* s ~ AX\ P - in) 



Vc b v: 



AX 3 , 



V ub V* s ~ AX 3 (p- in) 



(3.4) 



which changes the relative importance of the corresponding hadronic amplitude parameters 
for BRs and CP asymmetries in comparison with the U-spin analysis of D° decays. The 
latter can be parametrized analogously to eq. (2.7), 



A[B°- 
A[B° 

A[B°- 
A\B° 



D+K~] 
DJD+] 



2 VcbV* Bw= 



1 



%=i 
Bw=i 



(X u + X c ) re^ + ne 



101 



(X u + A c ) (re 



r\ e 



2 VibV* B w= 



+ (A U -A C ) (l + r e^° 
-(A u -A c ) fl-roe** 1 



(3.5) 



where for simplicity we have used the same notation for the amplitude ratios as for the U- 
spin analysis of D° decays. As for the experimental input, we consider the four observables 



4 A similar analysis could be performed for U-spin related modes, where the B^-meson in the initial state 
is replaced by B s , and/or the effective Hamiltonian for 6 — > sqq' is replaced by the one for b — > dqq' . 
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[49] 

, „ Br\B° ->■ D+K-}/\vdk\ 

bn/br 2 = — ^ J , -T a 11.8 ± 3.6 , 

11 Br[£?° -+ K-7r+]/\p K7r \ 

, „ BrLB -»■ D-D + ]/\p D D \ 

_ Br[fl°-> P^+l/lfc, | 
4/ Br[B» -t Jr-jr+]/|ftf»| 1 ' 

and 

aC p 2 = ^ [5° K"vr + ] ~ -0.098 ± 0.013 . (3.7) 

From this we can perform a fit to the 8 real parameters (r, rQ,ri,r[ and 0, (f>o,<f>i, 4>i), where 
again we generate random points which fulfill the experimental constraints within at least 
2a and lead to a total x 2 -value less than 6. To reproduce the difference between the two 
tree decays B° -> D+R- and B° -)• L>,7 7r + , in principle, one can allow for two different 
classes of solutions, with either > 1 or r± < 1. In the following, we will focus on the case 
r[ < 1, which would be the natural choice in the factorization approximation. 

Despite the fact that we might expect W-spin to be badly broken, the parametrization 
is sufficient to reproduce the available experimental data. A comment is in order about 
the fine-tuning in the amplitude parameters which is necessary to suppress the BR for 
the decay B° — > K~ir + with respect to B° — > D~ D + . To quantify this effect, we have 
determined the box-counting dimension (bed) of the fitted parameter space, following [53]. 
We found a bed which is indeed somewhat smaller than the naive dimension, 4 (refering 
to 8 fit parameters minus 4 experimental constraints). While each individual contributing 
amplitude is basically unconstrained in the fit, it turns out that the result for the average 



r 2 + rf + Tq 



is bounded from above and, at the same time, again sets the order of magnitude for the 
potential direct CP asymmetries. This is illustrated in fig. 5, where we show results for the 
direct CP asymmetry 

acp 3 = ^cp[^° -> D s D+ \ (not yet measured) , (3.8) 

and for the ratio acp3-acp2 as a f unc ti° n of r avg . In the approximation of universal strong 
phases in W = (W = 1) amplitudes, the latter would be given by ^^t^~ ~ — 1- The 
following observations can be made: 

• The points generated by the fit satisfy 

0.6<r avg <6.0 (forx 2 <6), 

0.8 < r avg < 5.3 (for X 2 < 1)- (3.9) 
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Figure 5. Direct CP asymmetries as a function of r avg (3.8) from a fit of W-spin parametrization 
to 4 observables in non-leptonic B° decays. Plot conventions as in fig. 1. 



• The maximal size of the (yet unmeasured) direct CP asymmetry |acp3| in B° — > 
DjD + is given by ~ 12% for x 2 < 6 (restricting ourselves to values of x 2 < 1 we 
would get slightly reduced values ~ 8 — 10%). Of course, smaller values for acp3, 
as expected in the factorization approximation, are not excluded in our approach. 
Experimental measurement of the CP asymmetries in B° — > D~ D + would be very 
useful in this respect. 

Of course, we should keep in mind that this picture could change when higher-order W- 
spin breaking effects are taken into account. Still, we find this somewhat academic exercise 
useful in getting an idea about an upper bound on the generic amplitude ratios in non- 
leptonic B° and D° decays. Apparently, even for a badly broken symmetry like W-spin, the 
amplitude ratios in the considered B° decays do not exceed a value of 5 — 6 or so. Naively 
translated to the D° case, this implies that the amplitude ratios needed to explain the 
central value of AA*p in the SM are at the upper limit of the acceptable (or expectable) 
range. 
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4 Conclusions 



Let us summarize our main results: 

• The ratios of BRs in D° — > 7r + 7r~, ir + K~, K + K~ require 0(1) U-spin violating ef- 
fects. The required magnitude of U-spin violation can be understood (but not un- 
ambiguously predicted) from long-distance strong-interaction effects. 

• The relative strong phase between the U-spin symmetric and U-spin violating contri- 
butions has to be large, too, which points towards essential non-perturbative hadronic 
rescattering effects. This is not surprising from experience on dealing with D-meson 
decays. 

• As a consequence of U-spin violation, the direct CP asymmetries in D° — > ir + -ir~ and 
D° — > K + K~ are no longer related to be equal in magnitude and opposite in sign 
(neither in the SM, nor in NP extensions). 

• Within the SM, the hadronic matrix elements of c — > uqq operators with highly 
Cabibbo-suppressed CKM factors should be enhanced compared to the leading oper- 
ators by a factor of (3 — 5) in order to yield the observed central values for the direct 
CP asymmetries. Although we lack a comprehensive dynamical model to generate 
such an enhancement, there is no good reason to exclude the possibility of such nu- 
merical factors (the chiral enhancement of certain penguin operators in non-leptonic 
S-meson decays is a well-known example; indeed, in K decays such enhancements 
are even more pronounced). As a toy example, we have also studied the breaking of 
W-spin (u o c) in non-leptonic B° decays, which is found to exhibit a very similar 
pattern of amplitude ratios and strong phases. 

• Due to the presence of U-spin violation, NP contributions to the measured Avl^p are 
possible with both, NP operators having U = or U = 1. 

• Considering the specific model of a SM extension by a fourth generation, we stress 
that large parametric enhancement of direct CP asymmetries (i.e. with less suppres- 
sion in terms of the Wolfenstein expansion) in charm decays are not possible as a 
consequence of the tight constraints on the 4G mixing angles and CP phases from 
kaon and S-meson observables. Still, the additional short-distance contributions of 
the 4G quarks to the weak effective Hamiltonian allow for a numerical enhancement 
(but, in principle, also to a reduction) compared to the SM. 

• From the experimental point of view, it should be worth looking into other non- 
leptonic D-meson decay modes which could be accessible to LHCb or Super-B facto- 
ries. On the one hand, there can be modes like D + — > (f)ir + , D s — > 4>K + which 
are induced by the same operators in the weak effective Hamiltonian as D° — > 
ir + ir~ , K + K~ , and therefore could be expected to yield direct CP asymmetries of 
similar magnitude. As an example, for D + — > 4>tt + , given the BR of about 3.1 • 10 -3 
(including the BR for the analyzing decay <j) — > K + K~), with O(10 9 ) D mesons, a 
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5cj reach for a 0.5% asymmetry is possible. On the other hand, one would like to 
constrain direct CP violation in tree-level decay modes such as D + — > K°(K*°)ir + , 
D s —> (j)ir + (K + ) etc. in order to test against NP contributions in charged flavour 
transitions. 

• Briefly, it may be worth pointing out that analogous D° decay modes into light 
vector and pseudoscalar mesons (PV) can be included, specifically K*^K T , p ±r K T . 
It is readily seen that significant U-spin violation takes place: phase-space corrected 
BRs (similar to obsi in eq. (2.5)) give 0.29 and 0.21, respectively. These modes can 
therefore be used to complement our understanding of the important issues related 
to direct CP violation in charm decays. 

In conclusion, we have emphasized that the breaking of U-spin symmetry between 
strange- and down-quarks points towards large non-factorizable long-distance effects in 
non-leptonic D° decays with large strong phases. As a consequence, a SM interpretation 
of the present data on direct CP asymmetries in D° — > P + P~~ is plausible. On the other 
hand, NP models with a constrained flavour sector, like models with a fourth fermion 
generation, can contribute with a similar magnitude, leading to a moderate enhancement 
(or reduction) compared to the SM. 
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